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ABSTRACT 

The redshift at which the universe was reionized is currently unknown. We 
Q\ | examine the optimal strategy for extracting this redshift, z Te { on , from the spectra of 

early sources. For a source located at a redshift z s beyond but close to reionization, 
(1 + Zreion) < (1 + z s ) < + z re i on ), the Gunn-Peterson trough splits into disjoint 
Lyman a, (3, and possibly higher Lyman series troughs, with some transmitted flux 
in between these troughs. We show that although the transmitted flux is suppressed 
considerably by the dense Lya forest after reionization, it is still detectable for 
sufficiently bright sources and can be used to infer the reionization redshift. The Next 
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, Generation Space Telescope will reach the spectroscopic sensitivity required for the 
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1. Introduction 

The standard Big Bang model predicts that the primeval plasma recombined and became 
predominantly neutral as the universe cooled below a temperature of several thousand degrees at 
a redshift z ~ 10 3 (Peebles 1968). Indeed, the recent detection of Cosmic Microwave Background 
(CMB) anisotropics rules out a fully ionized intergalactic medium (IGM) beyond z ~ 300 (Scott, 
Silk & White 1995). However, the lack of a Gunn-Peterson trough (GP, Gunn & Peterson 1965) 
in the spectra of high-redshift quasars (Schneider, Schmidt & Gunn 1991) and galaxies (Franx et 
al. 1997) implies that the IGM is highly ionized at low redshifts, z ^ 5. Taken together, these 
two observations indicate that the IGM was reionized in the redshift interval 5 ^ ^reion £ 300. 
The epoch of reionization marks the end of the "dark ages" during which the cosmic radiation 
background was dominated by the ever-fading CMB at all wavelengths. During this epoch, ionized 
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hydrogen (HII) started to occupy again most of the volume of the universe. Most likely, the phase 
transition from HI to HII was triggered by the emission from the first stars and quasars in the 
universe (see Loeb 1998, and references therein). 

The characteristic distance between the sources that ionized the universe was much smaller 
than the Hubble scale, and the overlap of the HII regions surrounding these sources was completed 
on a timescale much shorter than the Hubble time. The resulting reionization redshift, z ie i on , 
is therefore sharply defined and serves as an important milestone in the thermal history of the 
universe. Measurement of this parameter would also constrain the small-scale power of primordial 
density fluctuations that produced the first objects. Theoretical models and three-dimensional 
simulations predict that reionization occurs at 7 ^ z r cion ^ 20, just beyond the horizon of current 
observations (Haiman & Loeb 1997, 1998; Gnedin & Ostriker 1997). Forthcoming instruments, 
such as the Space Infrared Telescope Facility (SIRTF), and the Next Generation Space Telescope 
(NGST) are expected to reach the sensitivity required for the detection of sources beyond the 
reionization redshift. It is therefore timely to ask: how could one infer 2; rc ; on from the observed 
spectra of these sources? 

The spectrum of a source at a redshift z s > z Ic i on should show a GP trough due to 
absorption by the neutral IGM at wavelengths shorter than the local Lya resonance at the 
source, A G bs < A Q (1 + z s ). By itself, the detection of such a trough would not uniquely establish 
the fact that the source is located beyond z rc ion ; since the lack of any observed flux could be 
equally caused by (i) ionized regions with some residual neutral fraction, (ii) individual damped 
Lya absorbers, or (iii) line blanketing from lower column density Lya forest absorbers (see, e.g. 
the spectrum of a z = 5.34 galaxy taken by Dey et al. 1998). An alternative approach proposed 
by Miralda-Escude (1998) is to study the detailed shape of the damping wing of the GP trough. 
Although a measurement of this shape could ideally determine the optical depth of the absorbing 
HI slab between the source redshift and reionization, it is also likely to be compromised by 
contaminating effects such as: (i) redshift distortions due to peculiar velocities; (ii) uncertainties 
in the absorption profile due to the unknown intrinsic shape of the Lya emission line from the 
source; (iii) ionization of the IGM in the vicinity of the source due to the source emission; and (iv) 
the existence of a damped Lya absorber outside or inside the source. 

In this paper we consider a different approach for measuring z re ion- Our method relies on the 
existence of some transmitted flux between the Lya and Ly(3 absorption troughs in the spectra of 
sources which are located just beyond the reionization redshift. For such sources, the GP troughs 
due to different Lyman series resonances do not necessarily overlap, and the transmitted flux 
between these troughs is only affected by the post-reionization Lya forest. In §2, we simulate the 
absorption by the dense Lya forest (or "Lya jungle") and demonstrate that the transmitted flux 
features are detectable. In §3, we estimate the number of sources which are sufficiently bright to 
allow detection of these features with NGST. Finally, §4 summarizes the main conclusions of this 
work. 
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2. Spectra of High Redshift Sources 

The phase transition from HI to HII occurs almost simultaneously throughout the IGM 
volume, when the expanding HII regions around the separate ionizing sources overlap (Arons & 
Wingert 1972). This process is nearly instantaneous (Haiman & Loeb 1998; HL98) and defines 
a corresponding redshift, ^overlap- However, at ^overlap each individual HII region could still have 
a non-negligible Lya optical depth due to its residual HI fraction. The GP trough from the 
uniform IGM disappears only at a somewhat later redshift, z ve i on , when the average Lya optical 
depth of the smooth IGM drops below unity. (We distinguish between this GP absorption and the 
absorption produced due to IGM dumpiness by the post-reionization Lya forest). In this paper, 
we assume that the delay between .z over i a p and .Zreion is small and set z over -\ a p = .z re i on , based on the 
following reasoning. 

The Lya optical depth across a stationary HII sphere is as large as ~ 10 4 (Osterbrock 1974), 
but cosmological HII regions are much thinner because of the cosmological redshifting of the 
Lya photons away from resonance. If Overlap ~ 7 and the HII regions are driven by quasars or 
star-clusters in Cold Dark Matter halos of virial temperature of ~ 10 4 K (Mh a io ~ 1-5 x 1O 8 M0), 
then the typical proper radius of an HII region is -Rhii ~ 0.05 Mpc (HL98). When two such 
HII regions first overlap, the neutral fraction x = nm/nu at their intersection is x ~ 2 x 10 -5 , 
assuming photoionization equilibrium^. Taking this value as representative for the IGM, the 
resulting Lya optical depth of the smooth IGM is T a ~ 3. The GP trough disappears as soon as 
the neutral fraction is reduced by another factor of ~ 3, i.e. when the local UV flux is increased 
by the same factor. Such an increase occurs shortly after the overlap epoch, since it requires the 
cumulative contribution from sources in a region as small as ~ 3 x 0.05 Mpc = 0.15 Mpc, i.e. 
~ 0.1% the local Hubble length. Therefore, the resulting breakthrough redshift, z rc i on is only 
slightly lower than Overlap • The short redshift delay distorts the shape of the damping wings at 
wavelengths extending by a few tenths of a percent ( ^ tens of A) below the blue edge of the 
GP trough near the observed wavelength of A a (l + z re i on ). Our simplifying assumption of sudden 
ionization is therefore justified to this level of spectral resolution. Note that the Ly/3 optical depth 
is ~3.5 times below that of Lya, so that the Ly/3 and higher GP troughs are even more sharply 
defined. 

We now consider the spectrum of a source located at a redshift z s > z re ion- Lya absorption 
by neutral hydrogen between z re i on and z s suppresses the flux in the observed wavelength 
interval, (1 + -2 re ion)A a < A D b s < (1 + z s )X a ; Ly/3 absorption causes a second trough at 
(1 + -Zreion) A/3 < A ODS < (1 + z s )Xp, and higher Lyman series lines produce analogous troughs at 
shorter wavelengths. In the absence of other effects, the resulting spectrum would consist of a 
generic sequence of troughs separated by blocks of transmitted flux, as depicted in Figure |l|. In 
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this illustration we assumed that reionization occured suddenly as discussed above, at z rc i n = 7, 
but included the damping wings of the Lyman lines at both edges of the absorption troughs 
(see the Appendix in Miralda-Escude 1998). Note that a Lya photon emitted by the source can 
escape from its host's HII region and travel across the uniform IGM only if the HII region has 
an unusually large radius of at least ~1 Mpc. The emitted Lya photons are much more likely to 
be absorbed by the IGM. Below we consider observations on a telescope with a sufficiently high 
angular resolution, so that the diffuse Lya re-emission by the IGM can be ignored. 

The top panel of Figure || shows that the closer z s gets to Zreion, the narrower the GP troughs 
become and the more similar they look to individual Lya absorbers. Equation (|l]) below predicts a 
few damped Lya systems with neutral hydrogen column densities Nn *S 10 20 cm -2 in the spectra 
of a source at redshifts z s ^ 5. A single absorber with TVh = 10 20 cm -2 and b = 35 km s _1 would 
have a rest-frame equivalent width of 7. 7 A. In order to avoid confusion with such an absorber, 
the width of the GP troughs must be ^ 8(1 + z s )A. Note also that peculiar velocities could 
shift the location of the narrow troughs by ~ lOA. The bottom panel of Figure [l] also indicates 
that as the source redshift approaches the limiting value (1 + z s )/(l + z re ion) = A^/Aq = 32/27, 
only a single narrow block of transmitted flux remains between the Lya and Ly/3 troughs. At 
still higher redshifts this window disappears, leaving a continuous GP trough at wavelengths 
^obs < A a (l + Zreion) which carries no information about the reionization redshift. 

Reionization also leaves a distinct mark in the soft X-ray regime. The source spectrum would 
be strongly suppressed at the Lyman edge, A Q b s = 912(1 + z s )A, due to the large continuum 
optical depth in from H and He ionizations (r con t ~ 10 5 ). The optical depth drops at shorter 
wavelengths due to the sharp decline in the ionization cross-sections, a cx A 3 . The flux recovers 
its intrinsic amplitude at E ^> 0.1 keV, where the optical depth falls below unity. The observed 
photon energy where this occurs provides a measure of the reionization redshift, as long as the 
value of Q^h 2 is independently measured. Would AXAF be able to detect this spectral signature? 
For a source whose observed flux at E = 0.1 keV is 100 nJy, the AXAF count rate in the 0.08-0.15 
keV band would be 10~ 3 per secondP], sufficiently large to collect 100 photons in 10 5 seconds. The 
ability of AXAF to determine z re \ on relies on the existence of sufficiently bright sources beyond the 
reionization redshift. In the case where these sources are early quasars, their required infrared flux 
is ^ 10/uJy. As illustrated by Figure |I] below, such sources are expected to be very rare at z > 5. 

The spectra shown in Figure [l] are highly idealized in that they take account only of 
absorption by the homogeneous HI gas between z s and z re i on . In reality, the blocks of transmitted 
flux shown in Figure [l] will be suppressed due to absorption by the residual HI islands in the 
clumped IGM which produce the Lya forest after reionization. In order to simulate the effect 
of the dense Lya forest (or "Lya jungle"), we created mock catalogs of absorption lines with 
statistical properties constrained by observational data at z < 5. Press & Rybicki (1993) have 
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analyzed in detail the statistics of Lya absorbers observed between redshifts 1.5 < z < 4.3, and 
derived a convenient form for the mean number n(Nn,b, z)dN^dbdz of clouds with HI column 
density between Nr and iVn + dN^, Doppler velocity parameter between b and b + db, and redshift 
between z and z + dz, 

n(N R ,b, z)= 2.63 x W- U (l + Z y ( ) \{b) cm 2 km" 1 s, (1) 

\l(J i4 cm z J 

with 7 = 2.46, (3 = 1.43, and 

{b-b f 



p(b) oc exp 



26? 



(b > 0) (2) 



with J °° p(b)db = 1, and the best-fit values bo = 32 km s" 1 , b* = 23 km s -1 . We used 
a Monte Carlo approach to generate a catalog of absorbers along the line of sight to the 
source, with redshifts, column densities, and fr-parameters in the ranges < z\ < z re i on ; 
2 x 10 12 cm" 2 < jV H| i < 10 20 cm" 2 ; and b { > 0, for 1 < i < N ahs = J ndzdN u db. The probability 
distribution of these parameters was chosen according to the above equations. Although these 
equation are based on data from redshifts 1.5 < z < 4.3, in the absence of additional data at 
higher z, we extend their validity to all redshifts below z re ion • A more rigorous extrapolation 
of this equation to high redshift would be to use three-dimensional numerical simulations of the 
high-redshift IGM (e.g., Hernquist et al. 1996). However, even such a calculation would be highly 
uncertain due to the unknown evolution of the flux and spectrum of the UV background at z > 5. 

In computing the observed flux F from the source at a wavelength A Q b s we include absorption 
by the first nine Lyman series lines from each of the iV a b s (£ re ion) absorbers along the line of sight, 

F{\ ohs )=F J] n e -^t/( W i,n). (3) 
ae[l,9] ie[i,iv abs ] 

Here Fq is the continuum flux that would have been observed in the absence of Lya absorbers; 
T i,n = (7re 2 /m e c)A^H,iA n / n /6i is the integrated optical depth for the n th Lyman line of the i th 
absorber; and A n and / n are the rest wavelength and oscillator strength of this line. The profile of 
each line is given by the normalized Voigt function U(x, a), where x\ <n = (c/6i)[A bs/ (l+^i) — A n ]/A n ; 
a- 1D = r n A n /4-7r&i; and T n is the natural decay rate of the n th Lyman line (see Press & Rybicki 
1993 for more details). 

Figure |^ shows five examples of source spectra processed through the mock Lya jungle. As in 
Figure [j], we assume sudden reionization at z re i on = 7 but include the damping wings of the GP 
troughs. The several thousand absorbers along a typical line of sight yield a considerable mean 
optical depth, r ~ 4.5. Although the underlying continuum flux is not visible at any wavelength, 
the figure reveals a series of narrow features of transmitted flux - remnants of the processing by 
the underlying Lya jungle - to which we refer hereafter as "transmission features". The spectra 
contain many such features with central flux values that rise infrequently up to a significant 
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fraction (up to ~ 50%) of the original continuum flux. These features could be detected by an 
instrument with a sufficiently high sensitivity and spectral resolution. 

As an example, we show in Figure ^ a blow-up of the spectrum around the Ly a and f3 
GP troughs of a source with a redshift z s = 7.08 . The spectrum indeed contains numerous 
transmission features; these features are typically a few A wide, have a central intensity of a few 
percent of the underlying continuum, and are separated by ~ lOA. With a sensitivity reaching 1% 
of the continuum, the edges of the Ly/3 GP trough between 8201 and 8283 A could be identified to 
a ~ io A accuracy. A similar precision applies to the Lya GP trough, although the damping wing 
there suppresses the flux within ~ 10 A of the edges by a factor of ~ 3. The longest wavelength 
at which flux is detected towards the blue edge of the Ly/3 GP trough, at ~ 8201 A, provides the 
best direct measurement of the reionization redshift. In practice, this method only yields a lower 
limit to Zreiom since the first observed transmission feature could be offset from the true blue edge 
of the GP trough. The uncertainty in the measurement of A J g(l + z re ion) would be equal to the 
typical separation between the transmission features. Figure |3] shows that at a sensitivity reaching 
1% of the continuum, the separation would be small (~ 10A), leading to a fractional error of only 
~ 10/8200 ~ 0.1% in the measurement of z TC i on . Note that this accuracy is similar to the level 
allowed by peculiar velocities of a few hundred km/s, which could move either observed edge of 
the GP trough by ~ lOA. 

It is apparent from Figure ||| that the above procedure applies only as long as the typical 
separation between the transmission features is much smaller than the width of the GP trough, 
i.e. for (1 + z s ) ^ 1.01 x (1 + Zreion)- The bottom three panels in Figure |2| also show that the 
transmission features gradually disappear as the GP troughs of the Lya and Ly/3 resonances 
approach overlap at (1 + z s )/(l + z TC i on ) = X a /X/3 = 32/27. Hence, suitable sources must lie in the 
somewhat more restricted redshift range, 1.01 ^ (1 + z s )/(l + z TC i on ) ^ 1.17. Finally, we note that 
an extrapolation of equation (|l|) to even higher redshifts predicts that the mean optical depth of 
the Lya clouds increases strongly with redshift. For example, using our mock Lya cloud catalog, 
we find that at wavelengths just below A Q (1 + Zreion), the mean optical depth increases from 

T = 4.5 for Zreion = 

7 to r = 13.1 for 

^reion — 10. Similar numbers are obtained by extrapolating 
the expressions for the mean optical depth from Madau (1995). This implies that detection of 
the remnant transmission features is considerably more difficult at higher redshifts, since the 
number density of lines above a given sensitivity threshold drops as ~ exp(— r). In order to detect 
transmission features separated by roug hly ~ 10A, as in Figure 0, the spectroscopic sensitivity 
must reach the approximate fractions of ~ 10 _1 , 10 -2 , 10 -3 , 10~ 4 , and 10 -6 of the continuum flux 
for z rc ion = 6, 7, 8, 9, and 10, respectively. 

3. Abundance of Suitable High Redshift Sources 

The detection of a single bright source from the epoch, 1.01(1 + z re ion) < (1 + z s ) < 
1.17(1 + z re ion)) could in principle provide an unambiguous measurement of z re ion- How likely 
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is it to find such sources? The best candidates are high-redshift quasars, since in analogy with 
their z < 5 counterparts, they are expected to have a hard spectrum extending into the far UV. 
Alternative sources are Gamma-Ray Burst (GRB) afterglows, supernovae, or primeval galaxies. 
The advantage of GRB afterglows is that they are bright and have featureless power-law spectra; 
however, their abundance at high redshifts might be small. Supernovae and galaxies may be more 
abundant than either quasars or GRBs, but are likely to be faint and possess soft spectra in the 
relevant UV range. 

As the abundance of quasars at z > 5 is currently unknown, we must resort to an extrapolation 
of the observed luminosity function (LF) of quasars at z 5s 5. Such an extrapolation has been 
carried out by HL98 using a simple model based on the Press-Schechter (1974) formalism. In the 
model of HL98, the evolution of the quasar LF at faint luminosities and high redshifts is derived 
from three assumptions: (i) the formation of dark-matter halos follows the Press-Schechter theory, 
(ii) each dark halo forms a central black hole with the same efficiency, and (iii) the light-curve 
of all quasars, in Eddington units, is a universal function of time. As shown in HL98, these 
assumptions provide an excellent fit to the observed quasar LF at redshifts 2.6 < z < 4.5 for an 
exponentially decaying lightcurve with a time-constant of ~ 10 6 years. The model provides a 
simple and natural extrapolation of the quasar LF to high redshifts and faint luminosities. At the 
faint end of the LF, the number counts of quasars are expected to be reduced because of feedback 
due to photoionization heating that prevents gas from reaching the central regions of shallow 
potential wells. To be consistent with the lack of faint point -sources in the Hubble Deep Field, we 
imposed a lower limit of ~ 75 km s" 1 for the circular velocities of halos harboring central black 
holes (Haiman, Madau & Loeb 1998). 

Figure || shows the total number N(z, F) of quasars brighter than the minimum flux F, 
located within the redshift interval 1.01(1 + z) < (1 + z s ) < 1.17(1 + z), in a 16' x 16' field (16 times 
the field of view of NGST). How many suitable sources would NGST detect? For z rc i on = 7, the 
average transmitted flux shortward of the Lyman a and (3 GP troughs is a fraction exp(— r) ss 1% 
of the continuum. Around z ~ 7, we expect roughly a single ~ 10 3 nJy source per 16' x 16' field, 
based on Figure |. Such a source has an average flux of ~10 nJy blueward of A Q (1 + z s due to 
absorption by the Lya jungle. According to the method depicted on Figure ||, a determination 
of Zreion to ~ 1% accuracy would require spectroscopy with a resolution of R— 100. Requiring a 
signal-to-noise of S/N=3, the necessary integration time for a 10 nJy source with NGST would be 
roughly one hour^J This is a conservative estimate for the detection of flux, since it assumes that 
the flux is uniformly distributed across the spectrum, while in reality it would be concentrated 
in transmission features that cover narrow wavelength bins and yield a higher S/N ratio in these 
bins. We therefore conclude that NGST will be able to measure the reionization redshift to ~ 1% 
percent accuracy up to z Te i on = 7. The accuracy would degrade considerably with redshift for 



5 This result was obtained using the NGST Exposure Time Calculator, available at 
assumed a mirror diameter of 8m. 
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■2-reion •■ 

The total number of sources per 16' x 16' field, irrespective of their redshifts, is ~ 3 x 10 3 , 
~ 600, and ~ 90 for a minimum intrinsic flux of .Fo=10, 100, and 1000 nJy, respectively. In the 
most optimistic case, for which the reionization redshift is close to the peak of the relevant curve 
in Figure ^, up to ~ 25% of all detected sources in each image would lie in the redshift range 
1.01(1 + z) < (1 + z s ) < 1.17(1 + z), and reveal the features depicted in Figure However, since 
only a small fraction of all baryons need to be incorporated in collapsed objects in order to ionize 
the IGM, reionization is likely to occur before the redshift at which the source number-count 
peaks (zpeak = 3.5,4.5, or 6 at sensitivities of 10 3 ,10 2 , or 10 nJy). Pre-selection of sources in 
particular redshift bins could be achieved photometrically by bracketing the associated location of 
their Lya troughs with two narrow-band filters. 



4. Conclusions 

We have shown that the Gunn-Peterson effect breaks into individual Lyman series troughs 
in the spectra of sources located in the redshift interval 1 < (1 + z s )/(l + z re ion) < 32/27. 
Although the transmitted flux in between these troughs is heavily absorbed by the Lya jungle 
after reionization, the residual features it shows allow a determination of the reionization redshift, 
^reiom for sufficiently bright sources. For a single source, a reionization redshift z re \ on ~ 7 could 
be determined to a precision of ~ 1% with a spectroscopic sensitivity of ~ 1% of the continuum 
flux. A simple model for the abundance of high-redshift quasars predicts that a single source per 
16' x 16' field could allow this detection at the ~ 10 nJy spectroscopic sensitivity of NGST. It may 
also be feasible to probe reionization with ground-based telescopes. Using the Keck telescope, 
Dey et al. (1998) have recently detected a source at 2=5.34, and found no continuum emission 
below the rest-frame Lya wavelength. The authors were able to infer a lower limit of r = 1.2 on 
the optical depth of the Lya forest between 1050 and 1170 A to z=5.34, which falls just short 
of the 1.3 < r < 1.9 implied by equation (||) for this redshift and wavelength range. Therefore, 
a somewhat more sensitive spectrum of this object could potentially constrain reionization at 
z < 5.34. 

The signal-to-noise ratio in the determination of 2: re i on might be improved by co-adding the 
spectra of several sources. In the absence of Lya jungle absorption, the cumulative spectrum 
would approach the sawtooth template spectrum derived by Haiman, Rees, & Loeb (1998, Fig. 1) 
and would be isotropic across the sky. An alternative signature of reionization should appear at 
the soft X-ray region of the spectrum (E ~ 0.1 keV), where the continuum optical depth to H 
and He ionization drops below unity. AXAF could detect such a signature in the spectra of only 
exceptionally bright objects. 

In order to enable a measurement of 2: re i on at redshifts as low as 6, it is essential that the 
wavelength coverage of NGST would extend down to ~ 7\p = 0.7 jim.. Detection of the transmitted 
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flux features requires progressively higher sensitivities as z re i on increases because of the likely 
increase in the optical depth of the Lya jungle with redshift. The nominal integration time 
of about one hour for a ~ 10 nJy sensitivity (with R=W0 and S/N=3) would be sufficient to 
determine z re i on up to a redshift of ~ 7 with a precision of ~ 1%. Following the extrapolation 
of equation (|l]) to higher redshifts, the required sensitivity needs to be one or two orders of 
magnitude higher if z re i n ~ 8 or ~ 9. Note, however, that if Zreion ^ 10, then it will be more 
easily measurable through the damping of the CMB anisotropics on small angular scales (HL98). 
Near the planned spectroscopic sensitivity of NGST, most sources are expected to have redshifts 
lower than z re ion- However, suitable high-redshift sources could be pre-selected photometrically in 
analogy with the UV dropout technique (Steidel et al 1996), by searching for the expected drop in 
the flux at the blue edge (A a [l + z re ion]) of the Lya GP trough. 

This work was supported in part by the NASA ATP grants NAG5-3085, NAG5-7039, and the 
Harvard Milton fund. We thank George Rybicki for advice on the statistics of the Lya clouds, 
and Peter Stockman for useful information about NGST. 
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Fig. 1. — Spectra of five sources located beyond the reionization redshift (here taken to be 
-Sreion = 7), when all contaminating effects, such as peculiar velocities, the ionization of the IGM in 
the vicinity of the source, or absorption by the Lya forest, are ignored. 
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Fig. 2. — Spectra of the five sources in Figure |], when the absorption due to the post-reionization 
Lya forest is taken into account. 
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Fig. 3. — Blow-up of the spectrum of a source at z s = 7.08, assuming sudden reionization at a 
redshift z re i on = 7. The solid curves show the spectrum without absorption by the high-redshift 
Lya forest, and the dashed lines show the spectrum when the damping wings are also ignored. 
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Fig. 4. — Predicted number of quasars within the redshift interval 1.01(1 + z) < (1 + z s ) < 
1.17(1 + z), with a minimum intrinsic flux F , in a 16' x 16' field (sixteen 4' x 4' NGST fields of 
view). The observed quasar luminosity function was extrapolated to high redshifts using the model 
in HL98. 



